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Abstract 

We investigate the average entropy of a subsystem within a global unitary orbit of a given mixed 
bipartite state in the finite-dimensional space. Without working out the closed-form expression of such 
average entropy for the mixed state case, we provide an analytical lower bound for this average entropy. 
In deriving this analytical lower bound, we get some useful by-products of independent interest. We also 
apply these results to estimate average correlation along a global unitary orbit of a given mixed bipartite 
state. When the notion of von Neumann entropy is replaced by linear entropy, the similar problem can be 
considered also, and moreover the exact average linear entropy formula is derived for a subsystem over 
a global unitary orbit of a mixed bipartite state. 
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1 Introduction 

In 1978, Lubkin [15] proposed a method of approximating the average entropy for a subsystem of a finite¬ 
dimensional quantum system in a global pure bipartite state by expanding the entropy as a series in terms 
of the average traces of powers of the system's reduced density operator, but the convergence of this series 
was never established. However, the author of recent paper [6] found an exact closed form expression for 
the average traces, in which Dyer gave a characterization of the convergence of the series. 

In fact. Page conjectured in [18] that if a quantum system of Hilbert space dimension mn is in a random 
pure bipartite state, the average entropy of a subsystem of dimension m ^ n should be given by the simple 
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and elegant formula 


m — 1 
In ' 


( 1 . 1 ) 


5m,n — Hmn 


where := j is the fc-th harmonic number. The average entropy Sm,n in (1.1) is also served as a way 
of understanding the information in black hole radiation. This formula was first proved by Foong and 
Kanno [7] by using Fourier transform, and next by Sanchez-Ruiz [20] and by Sen [19] by using random 
matrix theory connected with generalized Laguerre polynomials. Some years later, Lachal [14] used a 
probabilistic approach to give a re-derivation of Page's formula. Recenfly, Zhang [23] has shown fhat if a 
quanfum system of Hilbert space dimension mn is in a random pure bipartite state, the average diagonal 
entropy of a subsystem of dimension m ^ n should be given by fhe simple and eleganf formula 


cD _ TT 

^m,n — ^mn 


Hn. 


( 1 . 2 ) 


Based on fhe above menfioned formulas, i.e. (1.1) and (1.2), he derives quickly fhat the average coherence 
of a subsystem of dimension m ^ n is given by S,^ „ — Sm,n = (see also [26] for another approach). 

We know that a random pure state can be generated by a unitary operator chosen uniformly according 
fo Haar measure }i. Thus the above problem can be equivalently described as follows: Consider a complex 
quanfum sysfem AB which consisfs of fwo subsystems A and B. For a given pure state Pab = \^ab) {^ab\> 

I S{TrB{UpABU^))dii{U) = {^a < ^ b ), 

where S(p) := — Tr(plnp) is the von Neumarm entropy. Along this line, in a very recent paper [3], 
Christandl et al computed exactly the eigenvalue distributions of reduced densify mafrices of multipartife 
pure stafe by employing symplectic geometric method. Here we ask: can one have an analogous formula 
for a given mixed bipartite state pnB? This question has been paid no or little attention to the best of our 
knowledge. In fhis note, we will make an attempt to determine the average entropy of a subsysfem along 
a global unifary orbif of a given mixed biparfite sfate. 

In fact, recently, many researchers studied various problems along a unitary orbit of a quantum state. 
For example, the total correlation attained between the subsystems of a bipartife quanfum sysfem is 
constrained if fhe bipartife system undergoes global unitary evolutions. The authors of Ref [11, 12] in- 
vestigafed related problems motivated by some considerations in the field of quantum thermodynamics. 
They have not only obtained the value of the maximal quantum mutual information (QMI), but also the 
maximum QMI state in the balanced bipartite quantum systems. Unlike the maximum QMI case, finding 
fhe minimum QMI sfate on the unitary orbit is more difficult than finding fhe maximum QMI stafe in 
general. Luckily, they completely solved this minimum QMI state for two-qubit case. Besides, Zhang and 
Fei investigated relative entropy and fidelity between two unitary orbits of fwo stafes, respectively [24]. 
They have obtained a lot of compacf expressions for some exfremal values under consideration. We have 
already known that any bipartite quantum state can be diagonalized under the global unitary conjugation 
but cannot be achieved in general under the local unitary conjugations. Because of this, Zhang et al consid¬ 
ered the fidelity between one bipartite quantum state and another one undergoing local unitary dynamics 
[25]. The problems are related to the geometric measure of entanglement and the distillability problem. 
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Besides the above mentioned works, there is also one, where Oszmaniec and Kus [17] estimated the frac¬ 
tion of noncorrelated states within a unitary orbit of a given mixed bipartite state. Based on the obtained 
result, they have proven that the fraction of noncorrelated states tends to zero exponentially fast with the 
dimension of the relevant Hilbert space whenever the purity exceeding some critical value. Consequently, 
a state within a global rmitary orbit of a given bipartife stafe is asymptotically a correlated one. Motivated 
by this and Page's average entropy formula, in the present paper, we will consider the calculation of av¬ 
erage entropy along a global unitary orbit of a given mixed bipartite state. On the one hand, our attempt 
made here can be seen as a generalization of Page's formula. On the other hand, the obtained result can be 
viewed as an estimate for the entanglement within a global imitary orbit of a given mixed bipartite state 
since a random state within a global unitary orbit of a given bipartite state is asymptotically a correlated 
one. In deriving our main results, we get a lot of by-products which is of independent interest. 

The paper is organized as follows. In Sect. 2, we introduce the notion of unitary orbit for quantum 
states. In calculating the average entropy of a subsystem along a global unitary orbit of a given mixed 
bipartife state, we need to calculate some integral (Lemma 2.1), which is very important in this paper, 
over unitary groups. Then, some implications are discussed in Sect. 3. Specifically, we have obtained 
the following resulfs: (i) For a bipartite quantum system, maximally mixed state can be represented by 
a uniform probability mixing of tensor products of two marginal states at each point within the global 
unitary orbit of any given mixed bipartite state; (ii) We use the relative entropy and fidelity as a figure 
of merit for correlation, and estimate the average correlation along a global unitary orbit from above 
and below since analytical calculation seems unavailable; (iii) A detailed research is given to the sum of 
average entropies of two subsystems within a global unitary orbit. We summarize the main contents of 
fhis paper in Sect. 4. Finally, we present the detailed proofs of Lemma 2.1, Theorem 2.2, Proposition 3.2, 
and Theorem 3.4 in Appendix. 


2 Average entropy of a subsystem along a global unitary orbit 


Define unitary orbit Up of a given quantum state p on a d-dimensional Hilbert space T-Ld as follows: 

Up-= {upU^ -.UeUiUd)}. 

Choose any e Up^g for a given bipartite state p^B arid = TrB(p(^g). Let T> = Trg, then O* = Cjlg. 
Assume that F := 0*0. That is F(X) = TrB(X) CiIb- 


In what follows, we compute the average entropy along the unitary orbit Up^g of Pab- By the definition 
of entropy, S(p^) = — Tr (p^ Inp^), which can be rewritten as: 


S(p:,) = - Tr (p^igln [u^T{UpabU^)u]) . 


Note that — Inx = 1(1 — x)”, and it follows that 


-In 


1 °° 1 / r IN' 

U^TiUpABU^)u\ = J2 -[^a^'^b-[u^T{UpabU^)u\) 

n=l ^ 


= ^ --LJ+ [f(LJTH'^) 


n—1 


u, 
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where T := (g) 1b/ dg — Pab with Tr (T) = — 1. Thus 

= E ixr [viUTU^)]” Udp{U)y 

where d = dAdg- Denote 

fl„ :=Tr (p J IJ+ T{UTU^) ” Udp{U) 


Hence 


We see from this formula 


r °° a 

In(d^) ^ / S{p'j^)dp{U) = E V 

n=l " 


lim a„ = 0. 

n—)-oo 


( 2 . 1 ) 


( 2 . 2 ) 


Theoretically, although the specific formula for can be obtained [4, 5, 22], however the specific form 
of will be rather complicated. We can truncate this series (via keeping the first n terms) to obtain an 
approximation about the average entropy. That is. 


In(d^) > J S{p'y^)dp{U) ^ fli + ifl 2 


1 1 

+ ■ ■ ■ H —Ufi 

n 


(2.3) 


The more larger n is, the more tighter lower bound is. But however, the computing about a„ is becoming 
very complicated when n ^ 3. As an illustration, we truncate first two terms as an estimate below. More 
challenging task is to determine the exact value of the above series. 


(i). For n = 1, we have [22] 


j U'^T{UTU^)Ud}i{U) 


dTr(r(lrf))-Tr(r) 


dTr(r)-Tr(r(lrf)), 


d(d2-l) ■Tr(X)ld+ E(d2-1)' 


Since T = 0*0, choose any orthonormal basis {|<py) ; j = 1,... fdg} for subsystem space of B, it follows 
that 

da 

r(Z) = E M.jZMt 

hj=l 

where M/y = 0 |<p,) {q>j\ = ML which implies that Tr (r(ld)) = dAd\,Tr (T) = d^dg- Thus 


Therefore, 


(ii). For n =2, 


J U^TiUTU^)Udp{U) = [(1 + ddB)!^ 0 1 b - (d + dgjp^g] . 

[(l+ddB)-(d + dB)Tr (p^^y 


fli = 


d2-l L 


da 


f IJ+ [F(ljrLf+)j Udp{U) = E / U''MijUTU\MjiMM)UTU^M,^Udp{U) 




E / U'^MijUTU^MjtUTU'^MuUdpiU). 
i,1,1=1-' 

In order to calculate above integral, we will need the following result, described as follows: 


(2.4) 


(2.5) 
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Lemma 2.1. It holds that 

j UAU^BUXU^CUDU^dpi{U) = m ■ + pi2 ■ BC + ^3 ■ CB + m ■ B + pis ■ C, (2.6) 

where the coefficients }ij{j = 1,... ,5) can be found in the Appendix A. 

Proof. See Appendix A. □ 


This leads to the following result: 

Theorem 2.2. The average von Neumann entropy of a subsystem within the global unitary orbit of a generic mixed 


bipartite state a density matrix on TLd = is bounded from below by a quantity: 

J S{p'ffidy{U) ^ ai + ■ ■ ■ + ^ (VneN), 


where 


a,, = Tr 


pju^ T{UTU^) ”!Jdfi(!J) 




In particular, the first two terms can be given specifically: 

— 1 r 


fli = 


d2-l L 


(1 + ddg) -{d + dg) Tr (p^gj 




-[g + h- 


Tr 


(pab) + ^ Tr (^PAsj , 


where f, g, h are given below: 


/ = 

^ = 
h = 


^ll^B ^1 (j ,J Tr (\ I ^)(^B _-r \2 

(rf2-l)(d2_4) (^A + dfiTr 2j+ 1) 

2d(d^-l)(d2 -1)(4-1) 

(d2-l)(d2-4) 

{d\-l){d\-A)dl 

(d2-l)(d2_4) • 

J S{p'ffidy{U)^-ln{l-a3). 


Moreover, 


(2.7) 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


Proof. See Appendix B. □ 

Remark 2.3. In particular, for dA = dg = 2, then d = 4, we have 

“> = I - i (do)' “2 = ^-5’'' (do) ■ (2 W) 

This implies that 

ln2> yS(p^)dp(!J) ^max|^-^Tr(p^B),ln5-ln(2 + 2Tr(p^g))|. (2.15) 

It is left open that the explicit computing of the following integral for a given mixed bipartite state Pab- 

J In (u^T{UpabU^)U^ dy{U). (2.16) 
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Remark 2.4. We have already known that for any super-operator S over 'Hd [22], 


/ 

iU(d) 


dfi{U)U^Z{UXU^)U = Tr (X) X. 


d{d^-l) 


d{d^-l) 


(2.17) 


Now let d = d/^dg and T-Lii = I-La ® T-Lb with dim(?f, 4 ) = dA and dim(?fB) = dg. Assume that X = pab / 
a density matrix on Ha Fixing an orthonormal basis {\^B,j) ■ j = ■ -/de} for Hb- Suppose that 

r(X) = TrB(X) (g) Ig. Then it can be rewritten as: 


r(X) = \^B,i) \^B,j) {^B,i\)- 


hi=^ 


Clearly r(ly^ 0 Ig) = dgl^ 0 Ig, implying that 

Tr (r(l ,4 0lg)) = d,44 and Tr (T) = {dAStjf = dj^ds- 

y=i 


From the above discussion, we see that 

dp{U)U^Y{UpABU^)U = ® 1b + ^^^^PAB- 

Denote p'j^g = UpabU^ and = Trg(p(^g). Then 

Tr((p:4)') = Tr((p:,0lg)p:,g) =Tr(r(p:,g)p:,g) 

= Tr (D+F(Dp,4B!i^)!JpAB) • 


(2.18) 


Therefore 

(Tr ((p:4)')) := / MU)Tr(^{p'^)^) =Tr (^J dp{U)U^T{UpABU^)UpAB^ ■ 

That is, 

(Tr((4)^))^M + ^Tr(,i,). ,2.19, 

If we use the notion of linear entropy Sg(p) := 1 — Tr (p^), then we have seen from Eq. (2.19) that the 
average linear entropy is given precisely by 

/ SBpU)MU) = '■‘'^-^^1“^-^'! + ^J^SBpab). (2.20) 

We see from this that 

/ SM)MU) = + '-%^SdPAB). (2.21) 

Moreover for a mixed state Pab> f Sl(Pa)^P(^) — I SL(Pg)dp(li) if and only if dA = dg- We also see that 
J {SUp'a) + SL(p'g))dp(D) = + ^|^Sl(p,4b)- (2-22) 

In particular, if dA = dg = 2, then 

J SMdpiU) = J SMdpiU) = i + IsiipAB)- (2.23) 
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3 Some implications 


Some consequences of Lemma 2.1 can be presented as follows: Let O = Trg and Y = Tr^^. Then <1>* = (gjlg 
and Y* = 1^0. Thus for 


Pa ® p'b = ^*^{UpABU^n*^{UpABU^) (3.1) 

= E i'^A<^\i){i\)UpABU'^{\p){v\<^\j){i\)UpABU^{\v){p\(^tB), (3.2) 

i,j—l 

and it follows that 


lipA<^PB)^PW = E E (1a® |z)(il) f / UpABUH\p){v\<E>\i){i\)UpABU^du) {\v){p\(g)lB 

•’ i,j=l }i,v=l / 


leading to the following identity 


|(p>p'g)dH(Ll) = E^^. 


(3.3) 


This amounts to say that 


Proposition 3.1. For a bipartite quantum system, maximally mixed state can be represented by a uniform probability 
mixing of tensor products of two marginal states at each point within the global unitary orbit of any mixed bipartite 
state. 


Besides, we also get that 

J U\p'A^p'B)Udp{U) 

= E E [ ^\'^A<»\i){j\)UpABU^{\p){y\<»\j){i\)UpABU^{\v){p\<E>tB)Udp{U), 

iy=\ fi,v=i ■' 


leading to the following identity 

J U^{p'A<»PB)^dii{U) = Co - Id + Cl ■ Pab + C2 ■ p\b' 

where 

(d2 -l)(d|-l) f 2 W 

° (d2-l)(i2_4) 2Tr(^p^Bjj, 

d^(d^ + dg — 6) + 4 
^ (d2-l)(d2-4) ' 

_ 2d(d2 -l)(d2-l) 

^ (d2-l)(d2_4) • 


(3.4) 


(3.5) 

(3.6) 

(3.7) 


The total correlation attained between the subsystems of a bipartite quantum system is constrained if 
the bipartite system undergoes global unitary evolutions. The authors of Ref [11, 12] investigated related 
problems from the field of quantum thermodynamics. Since knowing the maximal possible variation in 
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correlations is useful, it raises the optimization problem, where a search of the maximal and minimal 
correlated states on a unitary orbit is needed. This is completely solved for two-qubit systems. 

In the following, we make an attempt to calculate the average quantum mutual information (QMI) 
along a unitary orbit of a generic mixed bipartite state. We provide some analytical upper and/or lower 
bounds for the average QMI and quantum fidelity along a unitary orbit, although we carmot obtain precise 
formulas for them. 

Proposition 3.2. A lower bound for the maximal correlation (defined by relative entropy) within the global unitary 
orbit of a mixed bipartite state is given asfolloivs 

y J(A : B)p,dy{U) ^ S (pnel ko ■ Id + ci ■ Pab + ci ■ Pab) ' (3-8) 

where I(A : B)p := ^ ps), where S(p||t 7 ) := Tr (p(Inp — Inn)) is the relative entropy. Similarly, the 

lower and upper bounds for the average fidelity within the global unitary orbit of a mixed bipartite state Pab 
given as follows 

Co-FciTr (p3ib)+C 2Tr (p5ib) ^ J 'P{p'aB'P'a^ Pb)^^^) ^'P{Pab,co ■ Id + ci ■ pab + C2 ■ p\b)- (3-9) 

where F(p,£r) := Tr (^^y/pcTy/p) is the fidelity between two states p and cr. 

Proof. See Appendix C. □ 

Remark 3.3. Recently, in [9] the authors introduced a general measure of correlations for two-qubit states 
based on the classical mutual information between local observables. They focus on (classical) correlations 
between sets of local observables insfead of the quantum vs classical distinction. In this perspective, quan¬ 
tum states can be characterized as a whole by the average amormt of (classical) correlations between all 
pairs of local observables. Under some restrictions, the authors calculated the average mutual information, 
whose value depends on the state purity, and the symmetry of the correlations distribution. 

For the two-qubit case, we apply the spectral decomposition pab = UAU^ with A = diag(Ai, A 2 , A 3 , A 4 ), 
where 

(Ai,A 2 ,A 3 ,A 4 )'^ e A 3 = |(P1,P2/P3/P4)^ e : Pj ^ 0(V/ = 1,2,3,4),^ Py = ■ 

Due to the invariance under unitary conjugation, we have 



The lower bound in (3.9) becomes 

4 

Co -h Cl Tr (^A^j C 2 Tr (|a^^ = cq + X^(ciAy -h C 2 A?). 

We numerically check the bounds in (3.9) by choosing 10000 random points in the probability simplex 
A 3 and drawing the 3D scatter plot of ('P(A),co -l- ci Tr (A^) -|- C 2 Tr (A^) ,F(A,cq -|- ciA -|- C 2 A^)). For 
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clarity, its three-view drawings are also given by using the corresponding 2D scatter plots: 'P(A) vs. 
F(A, Co -h Cl A -h C 2 A^) in Figure 1(b); cq -F ci Tr (A^) -|- C2Tr (A^) vs. F(A, cq -F CiA -F C2A^) in Figure 
1(c); and 'P(A) vs. cq -F Ci Tr (A^) -F C2Tr (A^) in Figure 1(d). We also notice that the maximum of 
Co -F Cl Tr (A^) -F C2 Tr (A^) is even less than the minimum of F(A, co -F ciA -F C2A^). 



Figure 1: Two-qubif case, (a) 3D scaffer plof of ( P(A),co -F ci Tr (A^) -F C2Tr (A^) ,F(A,co -F CiA -F C2A^) 
). 2D scatter plots: (b) 'P(A) vs. F(A, co -F CiA -F C2A^); (c) Co -F Ci Tr (A^) -F C2Tr (A^) vs. F(A, co -F CiA -F 
C2A^); (d) V{A) vs. Co -F ci Tr (A^) -F C2Tr (A^). 


In fact, not only the value of fhe maximal QMI is obfained, buf also the maximum QMI state pmax is 
derived in [11, 12]. Specifically, Pab = UAU^ for = dg, where A = {Xj : j = 1,... ,d = d^ds} is fhe 
spectrum of Pab> 


max I{A : B)p, = I{A : B)p_ = ln(d) - S(A), (3.10) 

Pab^^pab 


where 


d 

pmax — ^j\^j) {^j 


(3.11) 


for any generalized Bell-sfate basis {| fly)}. Unlike the maximum QMI case, finding Pmin on fhe unifary 
orbit is more difficult than finding pmax in general. Since QMI varies on a imifary orbit completely depends 
on the sum of two marginal entries, we need to figure ouf fhe spectra of fwo reduced sfafes of all the 
states in the unitary orbit Up. The literature indicates that calculating this set of compafible reduced sfates 
with a given spectrum of a global bipartife seems unforeseeable. This is well-known "quanfum marginal 
problem", which is fully solved for fwo-qubif case [1] and theoretically for fwo-qubit case [13]. Specifically, 
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for two-qubit case, the solution of quantum marginal problem is given by the following: Mixed two-qubit 
state pab with spectrum A = {Ai ^ A 2 ^ A 3 ^ A 4 ^ 0} and the marginal states pA and ps exist if and 
only if minimal eigenvalues Aj^j^ and A^-^^ of the marginal states satisfy the inequalities [ 1 ]: 

fmin(A;^i^,A® J ^ A 3 + A 4 , 

j '^min + ^min > ^2 + ^3 + 2 A 4 , (3.12) 

[I'^min-'^mml ^ min (Aj - A 3 , A 2 - A 4 ) . 

With the help of this result, the value of the minimal QMI on a unitary orbit is derived for two-qubit case 

[ 11 , 12 ]: 


■ B)p^^ = h{M + A 2 ) + h{Ai + A 3 ) - S(A), (3.13) 

where Ai ^ A 2 ^ A 3 ^ A 4 ^ 0 and h{p) := p]np + [1 — p) ln(l — p) is the binary entropy function defined 
for p e [ 0 , 1 ], moreover two-qubit state pmin is given by 

Pmm= (3.14) 

i/=i 

Here \ij is a re-indexing of A^-(Ago = Ai, Agi = A 2 ,...) and {| i)}, {|;)} are qubit basis states for subsystems 
A and B. 

By Eq. (3.8) and the concavity of von Neumann entropy, we can draw the following conclusion: 
Theorem 3.4. For any generic mixed bipartite state pab, denoting 

{Sa + Sb) := J S{p'A)dp{U) + J S{p'B)dp{U), 

we have that 

^(Pab) + S (^Pab||co ■ Id + Cl ■ Pab + C2 ■ Pab) ^ (Ba + Sb) -H ci ■ pab + C 2 ■ Pab) ' (3.15) 

where Cj{i = 0,1,2) are given by (3.5), (3.6), and (3.7), respectively. Furthermore, 

S (|DAB||co'ld + ci ■ Pab + C2 ■ Pab) ^ J Kd-'- B)p,dp{U) ^ S (cg -Id + ci ■ Pab + C 2 ■ Pab) -^{Pab)- 
In particular, we get a universal entropy inequality: 

S (pubI ko ■ Id + Cl ■ Pab + C 2 ■ Pab) ^ S (cg ■ -h ci ■ pab + C 2 ■ Pab) - S(pab)- (3.16) 

Moreover, pab = cg ■ Id -t- Ci ■ Pab + C 2 ■ p\b only if Pab is maximally mixed state, i.e., Pab = ^ ® 

Proof. See Appendix D. □ 

This result tells us whenever the state Pab is not maximally mixed state, we can always find correlated 
states on a global unitary orbit of this state. Indeed, all quantum states are imitarily cormected to classical 
states, one-way or fully classical [16]. 
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In particular, for the balanced bipartite system, i.e., d/^ = dg, we have, for d = d/^dg, 

J^ipAWW + J S{pB)d}i{U) ^ max(S(p^)+S(pB)) =maxJ(A:B)p, + S(pAB)=ln(d), 

where we used the fact that max[j J(A : B)^/ = ln(d) — S(p, 4 b) [11]. Apparently, this upper bound is true 
for dj^ 7 ^ ds, and moreover it is trivially since S(p^) ^ ]n{dx), where X = A,B. Thus the following 
inequality is always true 

I S{p'A)dfi{U) + j S{p'B)dfi{U) ^ ln(d). (3.17) 

In view of this, we get a tighter upper bound for the sum of the average entropy of two subsystems: 

J HPa)^f{U) + J S{pB)dp{U) ^ S (co ■ Id + Cl ■ pab + C2 ■ Pab) ^ (3-18) 


We will check the difference between the upper bound S (cq ■ + Ci ■ Pab + ^2 ' Bab) lower 

boimd S(pab) + S (babI ko ■ Id + ■ Pab + ^2 ■ Bab) since we want to identify the range of the sum of two 

average entropy. Let 


•^(Bab) := S (^co ■ Id + Cl ■ Pab + C2 ■ Pab) “ S(bab) - S (babI ko ■ Id + ci ■ Bab + C2 ■ Pab) (3-19) 
and V = Tr (Bab)' which lies in [d~^, l] for d = dAdg- Apparently, J- ^ 0 over the whole set of states. 

If dA = ds = 2, then 

2-V 1 2 

' '‘ = 5- '^=5- 

Thus in the two-qubit case, we have 

J ld^(BA <^BB)^d/c(lJ) = ■ I4 -h ^ ■ Bab + g ■ Bab- ( 3 . 21 ) 


Note that S(pab) + S (bab I ko ' Id + ci ■ Bab + C2 ■ Bab) = - Tr (bab In [co ■ Id + ci ■ Bab + C2 ■ Bab] ) ■ Now 
J- is reduced to the following form: 


•^(Bab) 


2-V ,1 22 

-^■M + ^-Pab + ^-Pab 


+ Tr Pab In 


2-V 

10 


' + g ■ Pab + ^ ■ Pab 


(3.22) 


Since V is invariant under unitary conjugation, it follows, via the spectral decomposition p = UAU^ of 
BaB/ where A = diag(Ai, A 2 , A 3 , A 4 ), that 


JB(A) = S 


2-V 
10 ' 

Tr ( Ain 


l 4 + _.A+_.A^ 

2-V 


10 ■14+5-A+--A 

We see that is a symmetric function defined over the probability simplex 


^3 = \{pv V 2 , P3/ PiY e : Vj ^ 0(V; = 1,2,3,4), Py = 1 I . 


/=i 


(3.23) 
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Next we make numerical analysis of this function defined over the probability simplex. 


•^(^ 1 , A 2 ,A 3 , A 4 ) = 


^A^ln 

i=l 

'1-V 


l-V 1 


-k) 

(3.24) 

f2-P 1 2 

(3.25) 

10 +5'^^'+5^)' 




where (Ai, A 2 ,A 3 , A 4 )''' e A 3 and V = V{A) := Tr (A^) = Ey=i Ay. 

Choose random point A in the probability simplex A 3 , and then draw the 3-dimensional figure about 
the 3-tuples (x,i/,z) where x = 'P(A), y = — Tr ^Aln ■l 4 -l- 5 'A-|-§-A^ 'j, z = S(co -|- ciA -|- 

C2A ^). To illustrate the relationship among V{A), S{pab) + S {pabW^o ' Id + ■ Pab + ^2 ■ p\b) ^(^0 + 

Cl A -)- C 2 A^), we show the 3D scatter plot in Figure 2(a), where there are 5000 random points used in total. 

To make things more clear, we also demonstrate 2D scatter plots: S{pab) + S {pab I ko ■ Id + ■ Pab + ^2 ■ Pab) 

vs. S(co -I- CiA -h C 2 A^) in Figure 2(b), P vs. S(co -|- ciA -)- C 2 A^) in Figure 2(c), and P vs. P in Fig¬ 
ure 2(d). Note that for fhe case where Ay = |(/ = 1,2,3,4) we have S(A) = S(co-H CiA-|-C 2 A^) = 
ln4 = 1.39 and P = 0; while the case with one Ay equaling to 1 and the other three being zeros yields 
S(A) = 0, S(co-H CiA-h C2A2) = ^ Tr (Aln(co-h CiA-|-C 2 A^)) = hi^ and 

P = ~ 0-58. In these plots, P S [|, 1], S(A) e [0,ln4], S(co -h ciA -|- C 2 A^) S [0.94, In 4] 

and P e [0,0.58]. 



Figure 2: Two-qubit case (d ,4 = ds =2). (a) 3D scatter plot of (T’(A), — Tr (Aln(co-|-CiA-|-C 2 A^)) ,S(co-H 
ciA -|- C 2 A^) ). 2D scatter plots: (b) — Tr (Aln(co -F ciA + C 2 A^)) vs. S(co -F ciA -F C 2 A^); (c) P{A) vs. 
S(co-F CiA-F C2A2); (d) P{A) vs. P{A). 

Remark 3.5. Denote p' = UpU^ for a fixed state p. If U is such fhat I{A : B)pi > I{A : B)p, then we 
say that U enhances the correlation between A and B, otherwise, we say that it weakens the correlation 
A and B when I{A : B)pi < I{A : B)p. In particular, if 1{A : B)pi = 0, then we say that U decouples A 
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from B. Clearly there exists a state (for example completely mixed state) of which its correlations carmot 
be enhanced by any unitary. We see from previous discussion, finding minimal mutual information is in 
terms of an optimization problem that is extremely hard to handle in higher dimensions. Proposition 3.2 
and Theorem 3.4 can be seen as a different strategy based on probabilistic averages. Denote 

A{Pab) ■■= j I(A : B)p,dp(U) - I(A : B)p. 

If A(p) > 0, then the correlations existing in p, 4 B can be enhanced by many unitaries (via a perspective 
of concentration of measure phenomenon from Measure Theory). Clearly a lot of states of product forms 
satisfies A(p) > 0. In the opposite, A(p) < 0 means that the correlations existing in can be decreased 
by many unitaries. Naturally, some questions arise: For a given state Pab/ one wants to know wether if its 
unitary orbit of pab contains product state(s). 

Remark 3.6. A similar questions can be considered for the so-called quantum conditional mutual information 
(QCMl), defined by I{A : B\E)p := S{pae) + S(pbe) - S(pabe) - S(pe) for a tripartite Pabe on a tensorial 
Hilbert space T-La ® Td-B <8) 'He- We still denote = {U ^ l£)pABE(!d^ C) 1 e) for any unitary U on 
'Ha ®'Hb- Analogously, maximal and/or minimal QCMl can be considered, that is, 

maxt(A : B\E)pr and mint(A : B\E)pi. 

Clearly min(j I{A : B\E)pi is very important since it gives a lower bound for QCMl: 

I{A : B\E)p ^ rmnt(A : B\E)p,. 

We can also consider the average QCMl: 

J l{A:B\E)p,dp{U). (3.26) 

All problems mentioned above are beyond the goal of this paper. We will come back to them in the future 
research. 


Remark 3.7. Given a quantum charmel £ (trace-preserving and completely positive linear map) which is, 
via Kraus representation, represented as £ = Ey Adg^., where Ad£^(X) := EyXEt. We can use Lemma 2.1 
to get the average purity of a rmitary orbit of a given state undergoing a fixed quantum charmel £. Indeed, 
by Choi-Jamiolkowski isomorphism [2], /(E) := (E 0 l)(|ld)(ld|) = |Ey)(Ey|, 


J E(!Jp!J+)^dp(!i) 


^E,- (UpU^EjEjUpU^'j eJ 
cj 


d Tr (p^) — 1 
d(d 2 -l) 


ETr(ErE;) 



d — Tr (p^) 

d(d 2 - 1 ) 


E(ld)l 


(3.27) 

(3.28) 


implying that 

/Tr (£(UpU*)^) MU) (/(O^) + , (3.29) 

where the Ihs of the above last formula is just the average purity of a unitary orbit of a given state 
undergoing a fixed quantum charmel E. 
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4 Concluding remarks 


In this paper, we investigate the average entropy of a subsystem along a global unitary orbit of a given 
mixed bipartite state in the finite-dimensional space. Although it is still unable to derive the closed- 
form compacf formula for the mixed state case, compared with Page's formula, neverfheless we get an 
analytical lower bound for this average entropy for the mixed state case. In deriving this analytical lower 
bound, we obtain some useful by-producfs of independenf interest, for insfance, for a bipartife quantum 
system, maximally mixed state can be represented by a uniform probabilify mixing of fensor producfs of 
fwo marginal slates at each point within the global unitary orbit of any given mixed bipartife slate. This 
is amazing. In addition, from the discussion after finishing the proof of Proposition 3.2, we see fhat the 
average entropy of a subsystem is intimately related to the well-known "quantum marginal problem" or 
N-representability in quantum chemistry. Besides, it also connects with entanglement polytope, a notion 
proposed by M. Walter et al. in studying multipartite entanglement from single-particle information [21]. 

These obtained results can be applied to estimate average correlation along a global unitary orbit of 
a given mixed biparfite stale. The corresponding numerics about these results is also provided in lower 
dimensional case (for insfance d/j = dg = 2). Excepf thaf, fhe resulfs obfained in the present paper can 
also be used to study further the average coherence of a class of random stafes induced from isospecfral 
bipartife mixed stafes. Indeed, recenfly we have already calculafed exactly the average coherence for 
random mixed quantum states [23, 26] induced from random biparfite pure states. 

Finally, we conclude this section with two open problems: (i) computing the average coherence (via 
relative entropy of coherence) of a subsystem of isospecfral bipartite systems 

J {s{{p'A)di.s)-S{p'A))dp{U), (4.1) 

and (ii) proving the following identify for fhe balanced system, i.e., = dg, 

J S{p'A)dii{U) = J S{p'B)dp{U). (4.2) 

If (4.2) zvere frue, fhen we would see from Theorem 2.2 and Theorem 3.4 that 

- ln(l - fli) ^ J S{p'A)dp{U) < is (co ■ Id -6 Cl ■ pab + C 2 ■ p\b) ■ (4-3) 

We hope that the present work and questions proposed can bring out more interesting and insightful 
perspective(s) in quanfum information theory. 
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Appendix 


A. The proof of Lemma 2.1 

In this paper, we will utilize some notion of matrix integral [4, 5, 22]. The formula in Lemma 2.1 is given 
firstly. A detailed reasoning is presented here. 


The proof of Lemma 2.1. Firstly we note that 


h 


UAWBUXWCUDW 


'0 = E(^ 


h 


h'A 


AU^BUXW^CUD 






= E Ui^nU,j, {h\A\f^)(f2 


ii.n-nA.fvi’i 


4) v'l 


BUXWC 




ii.n-nA.fvi’i 


BUXWC 


iij {i2\U\j2) {h |D|ii) 

'2)02|D|/;) 




^hh 


X in I A| f 2 ) {i '2 |B| is) (is \U\ is) (is |X| f,) (i^ |Lf+| i^) (i^ |C| i 2 ) {h\D\j[) 


^hji ^hji ^hh ^Ai'i ^' 2/2 ^' 3/3 (ii 1^1 is) (is |B| is) (is |X|i^) (i^ |C| i 2 ) (i 2 |D|i(). 


h,h,h'h'hA>^'yfi'ifi5 


Then we have: 

J UAU^BUXW^CUDW^dU 


A; ;/B,v .• X,- ;/C, ,• D, 


hh4l'h,h,i'2,i'yi[,f2,A 


W : Lf,„„ LI,,,, W, ,v u.u. ,2 dU 


/ 1 J 2 h''3^h3j3 h '‘2 i 2 ,;l I / ^hlWl2]2^l3n^l{]Wl'2]Wiy3 


i2,i3,h'i2,h,i'2,i'yA,i2.A 


A; , 76,7 ,• X,' ,7 C,7 ,• D,- ,7 
7I/./2 72 //i 


X E (ill4(l))(i2l4(2))(i3l4(S))(/ll/v(l))(;2li^(2))(i3li^(3))Wg(cr7T 4 

\n,aeS2, j 

= E Wg(u7T-4 
n,a^S2, 


( 


\hhA42,hA'''yi'vi2’i3 


^h,i'2^i2'h^h,j3^A'^2^h,A ^^44(1)) (i2l4(2)) (i3l4(S)) (/ll4(l)^ (^214(2)) (/3|4(S)) 


In what follows, we compute this value step-by-step. 
(1). If {n,a) = ((l),(l)),then 


An,j^ifi3Xj„jfA,i2Dj2,A{h\A){W2) ihlQ ih\fi) ihlfi) (islis) 


i2,h,h'h,h,i2,i3'fvi'2'i3 


= Tr (AD) Tr(X)Tr(BC) (iili)). 


(4.4) 

(4.5) 
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(2) . If ( 7 r,c 7 ) = ((1), (12)), then 

^ (fi |f;) ikiQ ihlQ in Ifi) ihlfi) ihlfs) (4-6) 

h,h,h'h4i'i'2’iyivi'243 

= Tr {A) Tr (D) Tr (X) Tr (BC) {ii\i[). (4.7) 

(3) . If {71,0-) = ((1), (13)), then 

x: (fi If)) {h\Q {h\Q {j2\f2) {jslfi) (4.8) 

12.13, h42,h-i2-iyA-i2'f3 

= Tr (ADX) Tr (BC) {ii\i[). (4.9) 

(4) . If (71,(7) = ((1), (23)), then 

^ A^.3,^.,B,,, 3 X^ 3 ,^.,Q,(flit;) (z 2|4) (isks) ihlfi) HM) ihlfi) (4-10) 

12.13, hkSAi-A'i'ihh 

= Tr (DAX)Tr (BC) (zi|z;). (4.11) 

(5) . If ( 7 r,c 7 ) = ((1), (123)), then 

B,,, 3 X^ 3 ^^.,Q,, 30 ^ 3 ,^., (filfl) (( 2 |f^) ( 73 ^ 3 ) (/il; 2 ) ( 721 / 3 ) (73l7i) (4.12) 

12.13, h,h,i3,i'2h’A’f2’f3 

= Tr (A) Tr (DX) Tr (BC) ( 7 i|z(). (4.13) 

( 6 ) . If ( 7 r,c 7 ) = ((1), (132)), then 

Y, A^. 3 ,^.,B,,, 3 X^ 3 ,^.,Q,(filfl) (i2|7^) (73|73) ( 711 / 3 ) ( 72 I 71 ) ihlfi) (4-14) 

12.13, h,i2,i3,AA3,fl,i2'f3 

= Tr (D) Tr (AX) Tr (BC) ( 7 i|z(). (4.15) 

(7) . If (71,(7) = ((12), (1)), then 

Y A^. 3 ,^.,B,,, 3 X^ 3 ,^.,Q,(fill') (( 2 | 7 'i) ( 73 | 73 ) (/il7i) ihlfi) ( 731 / 3 ) (4-16) 

12.13.11.12.13, ^2,A'i'lhh 

= Tr(AD)Tr(X)(ii|BC| 7 l). (4.17) 

( 8 ) . If (71,(7) = ((12), (12)), then 

Y B,,,3X^3,^.,C,,,30^3,^., (fi 14 ) (( 214 ) (73I73) (/1I/2) (/2I71) (73I73) ( 4 . 18 ) 

k,i3,h,12,13,AAyfvi^Hl 

= Tr (A)Tr (D)Tr (X) {i, |BC| ((). (4.19) 

(9) . If ( 7 r,( 7 ) = ((12), (13)), then 

Y A^. 3 ,^., B,,, 3 X^ 3 ,^., Q,(fi| 4 ) (2214) (7314) (/ll/s) ihlfi) ihlfi) (4-20) 

12.13, h,i2,j3,i2'A'ii'i2'i3 

= Tr(ADX)(fi|BC|4). (4.21) 

(10) . If {n,a) = ((12),(23)),fhen 

Y A^. 3 ,^.,B,,, 3 X^ 3 ,^.,Q,(fi| 4 ) ( 2214 ) (7314) (/il 7 i) ( 72 I 73 ) ihlfi) (4-22) 

12.13, h,i2,j3,i2,iyi'vi2'i3 

= Tr(DAX)(2i|BC|4). (4.23) 
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(11) . If {n,o-) = ((12), (123)), then 

^ (fi 14) {i 2 \i[) ihir,) in \f 2 ) ihlfs) ihlfi) (4-24) 

h,h,h'h4?,'i'2’iyivi'243 

= Tr{A)Tr{DX){h\BC\i[). (4.25) 

(12) . If {n,a) = ((12), (132)), then 

E (fi 14) (2214) (^3|4) ik l/s) (klfl) (klfi) (4-26) 

12 . 13 , h'i2,h-i2-i3'fvf2'i3 

= Tr (D) Tr (AX) (z'l |BC| 4) . (4.27) 

(13) . If (7r,n) = ((13), (1)), then 

^ B,,,3X^3^^., Cy^Dy^., (2114) (2214) (^314) (h III) (;2|y2) (y3l;3) (4.28) 

12.13, h,h,i3,i'2,‘'3,ivi'2,i3 

= Tr (AD) Tr(X)( 2 i ICB| 4). (4.29) 

(14) . If (tt, (7) = ((13), (12)), then 

^ A^.3^^., B,y3X^3^^., Cy^Dy^., (2114) (2214) (^3l4) (H Ifl) ihlfi) iklfs) (4-30) 

12.13, h,i2,i3,i'2,i'3,fvi'2,i'3 

= Tr (A) Tr (D) Tr (X) (fj |CB| 4). (4.31) 

(15) . If (7r,n) = ((13), (13)), then 

Y, A^.3,^.,B,,,3X^3,^.,C,y3Dy., (2i|4) (2214) (^314) (Mis) ihlfi) (Dili) (4.32) 

12.13, h,i2,i3,i2,hil'hi3 

= Tr(ADX)(2i|CB|4). (4.33) 

(16) . If (7r,n) = ((13), (23)), then 

Y A^. 3 ,^., B,y 3 X^ 3 ,^., Qy^Dy, ( 21 14) ( 2214 ) (^3l4) iklfi) HM) iklfi) (4.34) 

12.13, h,i2,i3,i2,i'yi’i,i'2,i3 

= Tr(DAX)(2i|CB|4). (4.35) 

(17) . If (7r,n) = ((13), (123)), then 

Y A^.3,^., b,,,3X^3^^., (2114) (2214) (^ 314 ) in Ifi) ihlh) ihlf) ( 4 . 36 ) 

k,i3,h,i2,i3,i2,i'yfvi2'i3 

= Tr(A)Tr(DX)(2i|CB|4). (4.37) 

(18) . If (7r,n) = ((13), (132)), then 

Y A^.3,^., By3X^3^^., Cy^Dy, (2114) (2214) (^3|4) ik l/s) ihlfi) iklfi) (4.38) 

12.13, h,i2,i3,ik3'fvf2'f3 

= Tr (D) Tr (AX) {h |CB| 4) . (4.39) 

(19) . If (7T,27) = ((23), (1)), then 

Y By,3X^.3^^., (21 liO (2214) (^314) ill Ifi) ihlfi) iklf 3 ) (4.40) 

12.13, h,h,i3,i'2,i'3,ivi'2,i'3 

= Tr(AD)Tr(X)Tr(B)Tr(C) (2i|4). (4.41) 
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(20). If {n,o-) = ((23), (12)), then 


^'sl^) O’l Ifl) ihlfl) iklfs) 

h,h,h'h4?,'i'2’iyivi'243 

= Tr (A) Tr (D) Tr (X) Tr (B) Tr (C) {h\^^). 

(21) . If (7r,n) = ((23), (13)), then 

^ (filfl) {hV^) (f3|f^) (/il;^) (M/^) (/sill) 

i2h,h42,j3,i'2'‘'yii'i2'i3 

= Tr(ADX)Tr(B)Tr(C)(ii|i;). 

(22) . If (7r,n) = ((23), (23)), then 

(/llA) O^l/s) (/ 3 I/ 2 ) 

i2,h,h42,h,i2,iyi'vi2'i3 

= Tr(DAX)Tr(B)Tr(C)(/i|il). 

(23) . If (7r,n) = ((23), (123)), then 

^h’i2^i'2'^3^h’i3^iyh^h4'i *'^' 21 ^ 3 ) (^' 3 14) (/lli 2 ) O 2 I/ 3 ) (/3lil) 

i2,h,h,i2’h42,i'yi'i,i'2'j3 

= Tr(A)Tr(DX)Tr(B)Tr(C) (fi|4). 

(24) . If (7r,n) = ((23), (132)), then 

B,,, 3 X^ 3 ,^., Q,,, 30 ^ 3 ,(fi It;) (z2|4) {hVi) (Mis) (MA) {h\fi) 

i2,h,h'h4342,i'yj[,i'24'3 

= Tr (D) Tr (AX) Tr (B) Tr (C) (n |4). 

(25) . If (7r,n) = ((123), (1)), then 

Yj ^h4'2^i'2'h^h4'?Yiyi2^h'fi ^^ 11 ^ 2 ) ^'2|4) (/llA) O 2 I/ 2 ) Osl/s) 

i2,h,h’i243,i'2,iyi'i,j'2'i3 

= Tr(AD)Tr(X)Tr(C)(fi|B|4). 

(26) . If (7r,(7) = ((123), (12)), then 

Y ^h42^i'243^h43^i'342^h4'i ^^'44) (^’2|4) (^’314) (/il4) O 2 I 4 ) (/ 3 I/ 3 ) 

l2,l3,/l,/2,/3-'2-'3'/l'/2'/3 

= Tr(A)Tr(D)Tr(X)Tr(C)(fi|B|4). 

(27) . If (7r,n) = ((123), (13)), then 

^ B,,, 3 X^ 3 ,^., Q,,, 30 ^ 3 ,^., (fi|4) (2214) (^314) (Mis) O 2 I 4 ) (/3|4) 

h,l3./l//2,73,'2''3'/l'i2'i3 
= Tr(ADX)Tr(C)(2i|B|4). 

(28) . If (7r,n) = ((123), (23)), then 

Y ^h4'2^i'243^h4'3^i'342^i24[ 4'il4) (^214) (^3|4) (/ilA) 0'2|4) (/ 3 I/ 2 ) 

h,l3./l/;2,73,'2''3'/l'i2'i3 
= Tr(DAX)Tr(C)(2i|B|4). 


(4.42) 

(4.43) 

(4.44) 

(4.45) 

(4.46) 

(4.47) 

(4.48) 

(4.49) 

(4.50) 

(4.51) 

(4.52) 

(4.53) 

(4.54) 

(4.55) 

(4.56) 

(4.57) 

(4.58) 

(4.59) 
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(29) . If {n,o-) = ((123), (123)), then 

(/i 14 ) (z 2|4) ihlA) ihlfi) ( 72 I/ 3 ) (Mil) (4.60) 

h,h,h'h43’i'2’i'yivi'vi3 

= Tr (A) Tr (DX) Tr (C) {h |B| i[) . (4.61) 

(30) . If {K,a) = ((123), (132)), then 

^ A^.3,^.,B,, , 3 X^. 3 ,^., C,,,,3Dy, (fi|4) (i2|4) (i3l4) (iilis) (i2lii) (islii) (4-62) 

h,h,h,h,h,i2,iyi'vi'2'i3 

= Tr (D) Tr (AX) Tr (C) {h |B| i[) . (4.63) 

(31) . If (71,17) = ((132), (1)), then 

^ A^ 3 ,^.,B,,, 3 X^ 3 ,^.,Q,, 3 Dy, (fi|4) (2214) (7314) (Mil) (Mil) (Mis) (4-64) 

hhMl-hA’^'yi'vi'yk 

= Tr (AD) Tr (X) Tr (B) (ii |C| 4). (4.65) 

(32) . If (7r,i7) = ((132), (12)), then 

^ A^ 3 ^^.,B,,, 3 X^ 3 ,^.,C,,, 3 Dy, (fl|4) (7214) (^314) (illiM (i 2 |il) (Mis) (4.66) 

i2,h,h,h'h'i2’'yA’f2’f3 

= Tr (A) Tr (D) Tr (X) Tr (B) (fi |C| 4) . (4.67) 

(33) . If (7r,i7) = ((132), (13)), then 

Y, A^. 3 ,^.,B,,, 3 X^ 3 ,^.,C,,,, 3 Dy, (fi|4) (7214) (^314) (Mis) (Mil) (Mil) (4.68) 

i2,i3,h42,hAAyfi.i2'i3 

= Tr(ADX)Tr(B)(7i|C|4). (4.69) 

(34) . If (7r,i7) = ((132), (23)), then 

Y A^.3^^.,B,,,3X^3^^.,C,,,3Dy, (7 i| 4) (72|4) (73|4) (Mil) (i2li3) (Mil) (4.70) 

i2,h,h,h'hA2,iyj'l’f2’f3 

= Tr(DAX)Tr(B)(7 i|C|4). (4.71) 

(35) . If (7r,i7) = ((132), (123)), then 

Y A^. 3 ,^.,B,,, 3 X^ 3 ,^.,C,,,, 3 Dy, (7 i|4) (72|4) (73|4) (Mil) (Mis) (Mil) (4.72) 

hhAl^h.h.iyi'yfi.i'yA 

= Tr (A) Tr (DX) Tr (B) {h |C| 4). (4.73) 

(36) . If (7r,i7) = ((132), (132)), then 

Y A^. 3 ^^.,B,,, 3 X^ 3 ^^.,C,,, 3 Dy, (7i|4) (7214) (^314) (Mis) (i 2 lii) (Mil) (4.74) 

i2,h,h,12^3^2,i'yi'i.fyA 

= Tr (D) Tr (AX) Tr (B) {h |C| 4) . (4.75) 

Combing the 36 cases together gives the desired conclusion: 

J UAU^BUXU^CUDU*d}i{U) = ^ ■ + }i2 ■ BC + 113 ■ CB + pn ■ B + fis ■ C, (4.76) 
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where the coefficients }ij{i = 1,.. .,5) are given below: 


:= Wg(l,hl)Tr(AD)Tr(X)Tr(BC)+Wg(2,l)Tr(A)Tr(D)Tr (X)Tr(BC) 

+Wg(2,1) Tr (ADX) Tr (BC) + Wg(2,1) Tr (DAX) Tr (BC) 

+Wg(3) Tr (A) Tr (DX) Tr (BC) + Wg(3) Tr (D) Tr (AX) Tr (BC) 

+Wg(2,1) Tr (AD) Tr (X) Tr (B) Tr (C) + Wg(3) Tr (A) Tr (D) Tr (X) Tr (B) Tr (C) 

+Wg(3) Tr (ADX) Tr (B) Tr (C) + Wg(l, 1, 1) Tr (DAX) Tr (B) Tr (C) 

+Wg(2,1) Tr (A) Tr (DX) Tr (B) Tr (C) + Wg(2,1) Tr (D) Tr (AX) Tr (B) Tr (C), (4.77) 


^2 := Wg(2,l)Tr(AD)Tr(X)+Wg(l,l,l)Tr(A)Tr(D)Tr(X)+Wg(3)Tr(ADX) 

+Wg(3) Tr (DAX) + Wg(2,1) Tr (A) Tr (DX) + Wg(2,1) Tr (D) Tr (AX), (4.78) 

:= Wg(2,l)Tr(AD)Tr(X)+Wg(3)Tr(A)Tr(D)Tr(X)+Wg(l,l,l)Tr(ADX) 

+Wg(3) Tr (DAX) + Wg(2,1) Tr (A) Tr (DX) + Wg(2,1) Tr (D) Tr (AX), (4.79) 

Hi := Wg(3)Tr(AD)Tr(X)Tr(C)+Wg(2,l)Tr(A)Tr(D)Tr(X)Tr(C) 

+Wg(2,1) Tr (ADX) Tr (C) + Wg(2,1) Tr (DAX) Tr (C) 

+Wg(l, 1, 1) Tr (A) Tr (DX) Tr (C) + Wg(3) Tr (D) Tr (AX) Tr (C), (4.80) 

H 5 := Wg(3)Tr(AD)Tr(X)Tr(B)+Wg(2,l)Tr(A)Tr(D)Tr(X)Tr(B) 

+Wg(2,1) Tr (ADX) Tr (B) + Wg(2,1) Tr (DAX) Tr (B) 

+Wg(3) Tr (A) Tr (DX) Tr (B) + Wg(l, 1,1) Tr (D) Tr (AX) Tr (B). (4.81) 

We are done. □ 


Remark 4.1. Wg := Eai-Jc dto(^A) called the Weingarten function [4]. In particular, for A h 3, we 
have: 


Wg(l,l,l) = 

(f-2 

(d ^ 

(4.82) 

d(d2-i)(d2_4) 

1 d) N/ 

Wg(2,l) = 

1 


(4.83) 

(d2-l)(d2_4) 

Wg(3) = 

2 

2 1 

(4.84) 

d(d2-i)(d2_4) 

d N/ 
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where = {d^ — 1) {(fi — 4). With these coefficients, we then have 

Ndfii := (^d - Tr (AD) Tr (X) Tr (BC) Tr (DAX) Tr (B) Tr (C) + ^ Tr (A) Tr (DX) Tr (BC) 

+ ^ Tr (D) Tr (AX) Tr (BC) + ^ Tr (A) Tr (D) Tr (X) Tr (B) Tr (C) + ^ Tr (ADX) Tr (B) Tr (C) 

- Tr (A) Tr (D) Tr (X) Tr (BC) - Tr (ADX) Tr (BC) - Tr (DAX) Tr (BC) 

- Tr (AD) Tr (X) Tr (B) Tr (C) - Tr (A) Tr (DX) Tr (B) Tr (C) - Tr (D) Tr (AX) Tr (B) Tr (Cij4-85) 


[d - 


d 


Tr (A) Tr (D) Tr (X) + ^ Tr (ADX) + ^ Tr (DAX) 


- Tr (AD) Tr (X) - Tr (A) Tr (DX) - Tr (D) Tr (AX), 

Nd}i 3 := ^Tr(DAX) + ^Tr(A)Tr(D)Tr(X)+('d-^^Tr(ADX) 


Ndjii := 


NdFs ■= 


- Tr (AD) Tr (X) - Tr (A) Tr (DX) - Tr (D) Tr (AX), 

( d - ^ ) Tr (A) Tr (DX) Tr (C) + § Tr (AD) Tr (X) Tr (C) + § Tr (D) Tr (AX) Tr (C) 

- Tr (A) Tr (D) Tr (X) Tr (C) - Tr (ADX) Tr (C) - Tr (DAX) Tr (C), 

^ Tr (AD) Tr (X) Tr (B) + § Tr (A) Tr (DX) Tr (B) + ( d - ^ ) Tr (D) Tr (AX) Tr (B) 

- Tr (A) Tr (D) Tr (X) Tr (B) - Tr (ADX) Tr (B) - Tr (DAX) Tr (B). 


(4.86) 


(4.87) 


(4.88) 


(4.89) 


B. The proof of Theorem 2.2 

By Lemma 2.1, 

J U^MijUTU^MjiUTU^MnUd}i{U) = l)-td + g{i,i, 1) ■ T + h{i,i, 1) ■ T^, 

where 

/(f,;, 1) = Wg(l, 1, l)d\Sji Tr (t^) + Wg(2, l)d\S^jSjiSu Tr (t^) + Wg(2, Tr (t^) 

+Wg(2,l)d,iTr (t 2) +Wg(3)4d,yTr (t^) + Wg(3)4d,7 Tr (t^) 

+Wg(2, l)d\Sj, Tr {Tf + W^{3)dlSij5ji5n Tr {Tf + Wg{3)dASijSjidu Tr {Tf 
+Wg(l, 1, l)dATT {Tf + Wg(2, l)d\Sij Tr {Tf + Wg(2,1)4^,-; Tr {Tf , 
g{i,j,l) = Wg(3)4dy; Tr (T) + Wg(2, l)d\S,jSjiSu Tr (T) + Wg(2, Tr (T) 

+Wg(2,l)d,iTr (T) + Wg(l,l,l)4d,yTr (T) + Wg(3)4d,, Tr (T) 

+Wg{3)d\5j, Tr (T) + Wg(2, l)4^4y;<5,, Tr (T) + Wg(2, l)d,4 Tr (T) 
+Wg(2, l)dA Tr (T) + Wg(3)4d,y Tr (T) + Wg(l, 1,1)4^,, Tr (T), 
h{i,j, 1) = Wg(2, l)d\Sji + Wg(l, 1, l)d\Sij5ji5H + Wg{3)dASijSjiSii 

+Wg(3)d,4 + Wg(2, l)d\Sij + Wg(2, l)4d,7, 

+Wg(2, l)d\Sj, + Wg{3)d\S,jSjiSu + Wg(l, 1, l)dA<^47d,,- 
+Wg(3)dA + Wg(2, l)4<5y + Wg(2, 1)4<5,,. 
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Note that the meaning of the notation Wg(*) can be found in the Appendix. Thus for / := i,l), 

g ■= gi^'i' 0 / and h := we have 

/ = ([Wg(l,hl)+2Wg(3)]d2 + Wg(2,l)(d4 + d + d4))Tr(r2) 

+ (3Wg(2,1)^2 + Wg(3) {dd\ + d)+ Wg(l, 1, l)d4) Tr {jf, 
g = ([4Wg(3)+2Wg(l,l,l)]d' + 2Wg(2,l)(d4 + ‘i + rf4))Tr(T), 

h = 6Wg(2, l)d^ + [Wg(l, 1,1) + Wg(3)]d4 + [Wg(l, 1, 1) + Wg(3)]d + 2Wg(3)d4. 


Hence, for 7 = 1^ (gjlg/dB — Pab/ Tr (T) = — 1 and Tr (T^) = +Tr (4 b)' 


g 

h 


Therefore, 


implying that 


d(d2-4-4 + i) 

(d2-l)(d2_4) 


Tr 


(t^) 


d^(4 - 3) + 2(4 “ 4 + 1) y 

(d2-l)(d2_4) ' 


d^(4 -i)(4-i) / , „ / 2 \ 

fd2-llW2_4) [d-A+deTr 2J 

(dA-lf, 


+ 


(d2-l)(d2_4) 

(d^- 24 - 2 ) (4-1) 


(d2-l)(d2_4) 


2d(d2 - 4 - 4 + 1) _ 2 d{dA - 1)(4 - 1)(4 - 1) 

(d2-l)(d2_4) ’■I ! (d2_l)(d2_4) 

^^(4-5)+44 (4-1) (4-4)4 

(d2-l)(d2_4) (d2_l)(rf2_4) ■ 


J T{UTU^) ^Udii{U)= f-ld + g-T + h-T^ 


«2 = / + ^'Tr(p^Br)+/z-Tr 

= +«■ (i “ ) + '■ ■ (4 ^ i ) 

“ + *4 '’ 4 ) ■ 4 '' 4 ) ^ ^ 


(4.90) 

(4.91) 

(4.92) 

(4.93) 

(4.94) 

(4.95) 


(4.96) 


We make further analysis of the term a,, although we have already known the fact that lim^-^oo a„ = 0. Let 
(pp{X) = Tr (pX). Apparently, cpp is a positive unital linear mapping (in fact, it is a positive imital linear 
functional from the set of n x n Hermitian matrices to IR). It is easily seen that f{x) = x” is a convex 
function from ]R to IR. By using [10, Theorem 4.15, ppl47], we see that 


Un — 


> 


fpQ U^T{UTU^)U "dp{U)^ /(4r(!JT4)!J)dp(!i) 

J dp{U){cppof){U^T{UTU*)U) > J d}i{U)f[cpp{U^T{UTU^)U)] 
J dp{U) Tr (^p4T(!jr4)Ll) " ^ (^Tr U J dp{U)U^T{ 


iUTW)U 
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By (2.4), we have 

Tr (^p j dii{U)U^r{UTU^)U 

Thus 


dA-1 

d^-1 


(1 + dds) -{d + ds) Tr 


= ai. 


an ^ fl?. 


(4.97) 


Then 


where Sl{pab) e 


^ {dA-l){dB-l) ^ (rfA^ll)(rf±MsL(p^B), 


d + 1 


d^-1 


0 , 1 - 


. Clearly 


{dA - l){dB - 1] ^ ^ (iA - 1) (l + hi) ^ 


d + 1 


d + 1 


(4.98) 


Now 

f S{p'A)dp{U) = E^^Ef = - ln(l - « 1 ). (4.99) 

n=l n=l 

We see from the above lower bound, i.e., — ln(l — ai), that when the purity of Pab decreases, fli increases. 
Hence such lower bound will be tighter. 


C. The proof of Proposition 3.2 


Clearly, the first inequality is easily obtained 


J I(A : B)p,dp(U) ^ S (^PabII / In [u^ip'A^ Pb)^] dp{U)^ 

cq ■ Id + Cl ■ Pab + ci ■ Pab) • 


^ S 

= s 

In order to show the second inequality, note that, for any two density matrices p and cr, 

F{p,cr) ^ Tr {^/p^/a) ^ Tr (pa). 

Then, for p'^g = UpabW^, we have 

F(Pab^Pa ® p'b) ^ Tr [p'abP'a ® Pb) = Tr (pab^I^ (Pa ® Pb) Ff) , 


(4.100) 

(4.101) 


(4.102) 


(4.103) 


implying that 

j F(PAB^PA®PB)dp(!i) > Tr (pabU^ (Pa®Pb) Ff) dp(!J) = cq + ci Tr (p^ib) +C 2 Tr (pE) -(4.104) 
By the concavity of fidelity, we have 

J F(pE'PA®PB)dp(lT) = J F(pab, Lf^(pAOPB)Ff)dp(!J) ^ F(pab,co ■ + ci ■ pab + C2 ■ pis)-(4-105) 

Therefore 

Co + Cl Tr (p3ib) +C2 Tr(p5^g) < j V{p'^g,p'^(E) p's)d}i{U) < F(pab,co ' Id + H ■ Pab + C 2 ■ P^b). (4.106) 
This completes the proof. 
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D. The proof of Theorem 3.4 

Note that = UpAB^^- see from (3.8) that 

J 1{A : B)p,dp{U) ^ S (pabIIco -Id + ci ■ Pab + C 2 ■ pig) 


(4.107) 

Since 1{A : B)p, = S{p'^) + S(pg) — S(pab)/ it follows that 

J (S(pi)+S(pg)-S(pAB))dp(!J) ^s(pAB||C0'ld + Cl-pAB + C2-pig) • 

(4.108) 

That is, 

(Sa + Sg) ^ S(pab) + S (^Pab Cq ■ Id + Cl ■ Pab + C2 ■ pig^ 

)• 

(4.109) 

This confirms the first inequality. Besides, by Eq. (3.4), we get 

S (co - Id + Ci -pAB + C2 - pig) =^(^j (p'a^ Pb)^^pW'^ 

(4.110) 

^ Js(u\p'a^p's)u) dp{U) = J S{p'A)dp{U) + J S(p'g)dp(Lr). 

(4.111) 


This confirms the second inequality. Therefore we have 


S(|Dab) + S (|Dab||co -Id + Ci ■ Pab + C2 ■ p\b) ^ S [co ■ Id + Cl ■ Pab + C2 ■ p\b) ■ 

This is equivalent to the following 

s (pabI ko ■ Id + Cl ■ Pab + C 2 ■ Pab) < S (^co ■ Id + Cl ■ pab + C 2 ■ Pab) - S(pab)- 

Next we show that S (cq ■ Id + Ci ■ Pab + Ci ■ p\b) ~ ^(Pab) if arid only if p^B is maximally mixed state. 
Clearly if jOyiB is maximally mixed state, i.e., S(jDab) = lir(d), since S (cq ■ Id + Ci ■ Pab + C2 ■ Pab) “ ^{Pab) ^ 
0, then S (cq ■ Id + ci ■ pab + C2 ■ Pab) ^ Iri(d), apparently S (cq ■ Id + ci ■ pab + C2 ■ Pab) ^ thus 

S(Pab) = S (co ■ Id + Cl ■ Pab + C2 ■ Pab) ~ the maximum of von Neuman entropy. Reversely, if 

S (co ■ Id + Cl ■ Pab + C2 ■ Pab) ~ S(pab)/ then by the obtained inequality, we have 

S (pabIIco ■ Id + Cl ■ Pab + C2 ■ Pab) = 0, 

which holds if and only if pab = Cq ■ Id + Ci ■ Pab + C2 ■ p^ig- This means that for any eigenvalue A(^ 0) of 
Pab must satisfy that 

C2A^ + (ci — 1)A + Cq = 0. 

Solve this equation, we get 

(1-Ci)- a/(1-Ci)2-4coC2 1 

=- 2F2 -= T 

Note that we have dropped another root being larger than one. Thus Pab is maximally mixed state. In 
fact, we get that Pab = cq ■ Id + Ci ■ Pab + C 2 ■ P^ig if arid only if pab is maximally mixed state. 
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